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ABSTRACT 

Given a flurry of recent claims for systematic variations in the stellar initial mass function (IMF), 
we carry out a novel inventory of the observational evidence using different approaches. This includes 
literature results, as well as our own new findings from combined stellar-populations synthesis (SPS) 
and Jeans dynamical analyses of data on ~ 4500 early-type galaxies from the SPIDER project. We 
focus on the mass-to- light ratio mismatch relative to the Milky Way IMF, Simf, correlated against 
the central stellar velocity dispersion, er*. For the SPIDER sample, we find a strong correlation 
between <5imf and av, independently of assumptions on the dark matter (DM) profile. The overall 
normalization of <5imf, and the detailed DM profile, are less certain, but the data are consistent with 
standard cold-DM halos, and a central DM fraction that is roughly constant with a+. For a variety 
of related studies in the literature, using SPS, dynamics, and gravitational lensing, similar results are 
found. Studies based solely on spectroscopic line diagnostics agree on a Salpeter-like IMF at high 
cr*, but differ at low a*. Overall, we find that multiple independent lines of evidence appear to be 
converging on a systematic variation in the IMF, such that high-cr* early-type galaxies have an excess 
of low-mass stars relative to spirals and low-cr* early-types. Robust verification of super-Salpeter IMFs 
in the highest-er* galaxies will require additional scrutiny of scatter and systematic uncertainties. The 
combined implications for the distribution of DM are still inconclusive. 

Subject headings: galaxies: evolution — galaxies: general — galaxies: elliptical and lenticular, cD 



1. INTRODUCTION 

The stellar Initial Mass Function (IMF) is fundamen- 
tally important to understanding both stellar popula- 
tions and galaxies. The Milky Way (MW) IMF was orig- 
inally characterized as a powe r-law mass-dis tribution, 
dN/dM oc M- a , with a ~ 2.35 (|Salneterlll955h . and sub- 
sequently refined to flatten at lower masses (M <^ 0.5M©; 
IKroupdl200ll IChabrierll200l . 

Whether or not the MW IMF describes stellar pop- 
ulations elsewhere in the Universe cannot yet be said 
through direct star counts. There have been some in- 
direct observational hints of IMF variations, and ample 
theoretical motivation for the se, but no broadly co nvinc- 
ing evidence has emcrg ed fcf. lBastian et al.ll2010fl . 

This situation has recently changed, with a flurry 
of studies of early-type galaxies (ETGs; ellipticals and 
lenticulars) turning up indirect evidence for systematic 
IMF variations. These studies use models of stellar 
population synthesis (SPS) to fit integrated-light data 
(broad-band colors and spectroscopic features), and fall 
into two broad categories: "pure" SPS, and "hybrid" 
SPS+gravitating mass analyses. 

The pure analyses rely on spectral lines that are dif- 
ferentially sensitive to giant or dwarf stars. These in- 
clude the TiO feature at 6130 A, the Na I doublet near 
8190 A, the Ca II triplet near 8600 A, the Wing-Ford 
[Fe/ Hl band at 9915 A, and the Ca I line at 10345 A 
(e.g.. iCenarro et al.|l2003l: Ivan Dokkum fc Conrovl 120101: 
iSpiniello et al.ll2012t iConroy fc van Dokkumll2012l h ere- 
after CvD12: lFerreras et al.ll2012t iSmith et al.ll2012D . 

The hybrid analyses assume an IMF and infer a stellar 
mass-to-light ratio Y* using a more conventional SPS 



approach based on colors, and age- and metallicity- 
sensitive spectral lines. Estimates of T* are also derived 
using dynamics or gravitational lensing. Comparison 
of the independent results th e n yields a revise d IMF 



eg 



iRenzini fc Ciottil fl99l iCappellari et all 
iFerreras et all 120081: iTortora et all 120091 



201 

Grillo fc Gobatll2010HTreu et al.ll2010t iNapolitano et al 
2010Ll2011l:IAuger et al.ll2010tlGrilldl2010HThomas et al 



Spin iello et all 2011 : IDutton et al.l ~l2012b.a: 
Sonnenfeld et al.ll2012t IWeener et al.ll2012ft . 



One may characterize a revised IMF through its T* 
relative to a MW-disk IMF, 5imf = Y*/T*,]yiw (the 
"mismatch parameter"), where for reference we adopt 
the Chabrier IMF. Remarkably, almost all the above 
studies found "heavy" IMFs (^imf 3> 1) for the most 
massive ETGs. Less massive ETGs, and spiral galaxies, 
appear to have "nor mal/light" IMFs (<5tmf 1; e.g., 
i Bell fc de Jonal [2001 iBershadv et all l201lf ISuvu et alJ 
120121: IBrewer et al.ll2012A~ 

These IMF findings are both potentially revolutionary, 
and highly controversial, and demand further investiga- 
tion. In particular, with hybrid analyses there are lin- 
gering questions about degeneracies associated with the 
distribution of nonbaryonic dark matter (DM). The time 
is also ripe to inventory the results to date, and see if 
the apparent emerging consensus holds up under quan- 
titative, systematic comparison - which could provide 
pressing motivation for understanding the physical ori- 
gins of the tre nds. Comparisons were made for some hy- 
brid analyses ((Thomas et al.l 1201 it IDutton et al.l l20l2at 
iWegner et al.ll2012f ) . but nor for the pure SPS work. 

In this Letter we carry out such an inventory, while pre- 
senting our own novel results for a large sample of ETGs 
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for reference, following the dynamical+SPS analyses of 
iTortora et all (|2012t hereafter T+12). We focus on the 
trends in Simf with central stellar velocity dispersion, a*. 
o\t is widely considered as crucially connected to galaxy 
evolution, and unlike T*, is relatively independent of the 
bandpass, and of Simf- 

We describe our data and analysis methods in Sec- 
tion [2j We present our results in Section [3] and make 
literature comparisons in Section [4] We summarize the 
conclusions and outlook in Section [5] 

2. DATA AND ANALYSIS METHODS 

We apply a combination of SPS and stellar dynamical 
models to a sample of ~ 4500 giant ETGs, in the red- 
shift range of z = 0.05-0.1, drawn from the Spheroid's 
Panchromatic Investigatio n in Different Environm ental 
Regimes project (SPIDER; |L"a Barbera et~aHl2010D . The 
details of sample selection and methods can be found in 
T+12, and here we provide only an overview of our two 
T*-estimation methods. 

Our data include optical+near-infrared photometry 
[from the Sloan Digital Sky Survey (SDSS) and the 
UKIRT Infrared Deep Sky Survey-Large Area Sur- 
ve yfl high-quality measurements of galactic struc- 
tural parameters (effective radius R c and Sersic in- 
dex n), and SDSS central-aperture velocity dispersions 
ffAn- The SPS-based T+ values were derived by fitting 
iBruzual fc Chariot! (|2003t l models to the multi-band pho- 
tometry, assuming a Chabrier IMF. 

Our dynamical-mass estimates use spherical isotropic 
Jeans equations fitted to the OAp data. We also extend 
the T+12 analysis using two-component mass models: 
a Sersic-based stellar distribution following the if-band 
light, and a standard DM profile. The data do not con- 
strain both components simultaneously, so instead we 
ado pt a series of plausibl e assumptions about the DM 
(cf. ICappellari efalll2012h . 

Our DM models are based on the iNavarro et al.l 
(|1996l hereafter NFW) profile, with an adjustable de- 
gree of baryon-induced adiabatic contraction (AC). For 
the virial mass and concentration (M v ; r , C yj r ), we adopt 
mean trends for a WMAP5 cosmology (Macci o~et al.l 
20081). We ta k e an M v i r -stellar-mass (M*) relation from 
Moster et afl (|2010t M+10 hereafter). Each galaxy's 
mass mod el then has one free parameter, Y*, plus op- 
tional AC (jGnedin et alJl2004L G+04 hereafter), provid- 
ing our no-AC-NFW and AC-NFW base-models. 

We explore the sensitivity of our results to these 
assumptions by doing the analyses with the follow- 
ing alternatives: a) AC recipes of vary ing strengths 
(jBlumenthal et al.lll986HAbadi et al1l2010l ): b) M vir with 
a fixed value for the entire sample: M vir = 10 12 M Q , 
10 13 M ^, or 10 14 M^; c) W MAP3-based c vir -M vir re- 
lation |Macci6etaL| |2QQH); d) no DM is present; c) 
mild kinematic anisotropy, with (3 = —0.2 or +0.2; f) 
c v i T -M v i T relation altered to mimic a warm dark matter 
(WDM) cosmology, a ssuming different particle masses 
|Schneider et al.ll2012l S+12 hereafter). 

In order to study the mean trends of with veloc- 
ity dispersion, we construct "average" galaxies from the 
data. Dividing our sample into 10 (TAp-bins, for which 
we compute median values of M*, R c n, i?Ap, and oa p . 

1 http://www.sdss.org, http://www.ukidss.org 
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Figure 1. Galaxy sample properties, binned by velocity disper- 
sion: physical aperture radius, SPS-based stellar mass (Chabrier 
IMF), effective radius, and Sersic index. The points and error-bars 
show the medians and 25—75 percentile scatter. 

We show these values, and their 25-75 percentile scatter, 
in Fig. [U 

Our final analysis-products, for each galaxy bin and 
mass model, will be the SPS-determined T*,mw, the 
dynamically-determined T*, the inferred <5imFj and the 
inferred central DM- fraction, /dm = 1 — T*/Td yn - For 
homogeneity, we convert oa p t o a e (the value at R e ) , 
using the best-fitting relation in ICappellari et al.l (|2006f ) 
(also done for the literature results later). 



3. RESULTS: IMF AND DM FRACTION TRENDS 

Our main IMF results from SPIDER are shown in Fig- 
ureO The two thick black curves correspond to our stan- 
dard no-AC-NFW (solid) and AC-NFW (long-dashed) 
models, and the suite of alternative models are also plot- 
ted, as labeled. 

It is clear that the overall normalization of Simf is de- 
generate with the adopted halo model, as DM can be 
traded against stellar mass. However, for a given flavor 
of halo model, there is always a strong correlation be- 
tween Simf and cr* (the Simf-M+ trend is weaker). We 
quantify this correlation for each model with a log-log fit, 
reporting the best-fit parameters in Table [1] the typical 
relation is Simf 9S °e- 

Othe r models of po tential interest are low-density DM 
cores dBurkertl Il995h . and alternatives to DM (e.g., 
iMilgroml |1983[ ). These would imply similar <5iMF-cr e 
slopes to our limiting no-DM model. 

One could alter these slopes with additional model- 
tuning; e.g., with strong- AC at high-cr e and halo ex- 
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Figure 2. IMF mismatch-parameter, <5jmf = T^/T* MWi vs. velocity dispersion <r e , for the SPIDER sample. Top: results for several 
different mass models. Horizontal lines show reference values for Salpeter, Kroupa, Chabricr IMFs (top-to-bottom). Bottom: residuals 
relative to the fiducial no-AC-NFW result, for an expanded scries of mass models (see legend at right, explanation in the main text, and 
reference list for abbreviations). The overall normalization of 5imf depends on the model assumptions, but a steep relation between <5imf 
and <T emerges as a robust result. 



Table 1 

Best-fit parameters for the relation log <5imf = a, - 
for model suite (see main text). 
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pansion at low-a e , one could completely flatten out the 
trend. We currently have no a priori motivation for such 
variations, but they would be a way to preserve IMF 
universality. 

We next examining the overall IMF normalization, 
with the no- AC and AC-NFW cases bracketing the most 
plausible range of models. For reference, we show 6\mf 



predictions for several standard IMFs (Salpeter, Kroupa, 
Chabrier). We note that at a fixed IMF, age, and metal- 
licity, the T* and <5imf values are uniquely predicted, 
but the reverse is not true. A given <5imf result can im- 
ply multiple IMF solutions, particularly if one allows for 
mass-functions more complicated than a pure power-law. 

If we adopt a no-AC-NFW model, a MW-like IMF is 
implied for low-er c galaxies, and a Salpeter IMF for high- 
er^ For AC-NFW, the low-er c galaxies have sub-MW 
IMFs, and the high-<7 ones have IMFs intermediate to 
Kroupa and Salpeter. 

In all cases, extremely bottom-heavy IMFs (a ;> 2.6) 
are ruled out, on average. Even if one assumed no DM, 
such IMFs would violate the dynamical constraints on 
the overall mass-to-light ratio. 

Although resolving the remaining IMF-DM degener- 
acy will require more extensive analysis, we carry out 
a simple exe rc ise to provide initial clues, inspired by 
IDutton et al.1 ()2012bf ). As in that paper, we select 
the galaxies with mean central stellar surface-densities 
> 2500 Mq pc~ 2 , and analyze them the same as the 
full sample. The rationale is that such galaxies are the 
most star-dominated, and the least sensitive to DM un- 
certainties. The results are shown by the gray curves in 
Figure [H where as expected, the model curves are closer 
together. The implied <5imf normalization is fairly low - 
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Figure 3. Inferred DM fraction within I R c , vs. <r c , for the 
SPIDER sample with different modeling assumptions. The solid 
and dashed black curves are our default no-AC-NFW and AC- 
NFW models. The gray curves are for fixed Chabrier IMF, with 
isothermal and constant-T^yn models, bracketing the range of 
non-parametric models (solid and dotted curves, respectively; see 
T+12). With standard ACDM halos and a variable IMF, / DM is 
constant or mildly dec reasing with <r e ; wit h a fixed IMF, /dm is 
strongly increasing (cf. Thomas ct al. 2011 figure 16). 



similar to the full-sample AC results. We may then apply 
this IMF result to the full sample if we as sume no ad- 
dition al systematic <5imf~£* correlation (cf. iSchulz et al.1 
l2010h . 

Although this paper is primarily concerned with the 
IMF, we briefly examine some implications for the cen- 
tral DM content. Our main mass models are, by con- 
struction, fully consistent with current expectations for 
ACDM halo profiles, while also agreeing with the obser- 
vations for a plausible IMF range (<$imf ~ 0.5-2.0). 

We show the implied /dm within 1 R c in Figure El 
We find fairly universal values of /dm ^0.2 and ~ 0.5 
for the no-AC-NFW and AC-NFW models, respectively. 
These results are not altered appreciably in the alterna- 
tive models explored above. Note also that the high-E* 
test above prefers the AC-NFW model. 

The Figure also shows that if we adopted a constant 
IMF, then we would infer a strong increase of /dm with 
cr e . Such behavior has been invoked as the driver for the 
"tilt" of the ETC fundamental plane (see iTortora et al.l 
2009 and references therein) , but now with renewed com- 
parison to realistic DM halo models, we find that the tilt 
is driven in part by the IMF. Put differently, the ob- 
served Tdy n _ Ce relation is too steep to explain through 
standard DM models, and requires an additional factor. 
Again, the IMF variation is in principle degenerate with 
systematic variations in AC, which we cannot completely 
exclude. 

4. LITERATURE COMPARISONS 

We now inventory other results from the literature, and 
see if they agree with our SPIDER-based results, showing 
the comparisons in various panels of Figure |4j We begin 
with those studies that used similar hybrid approaches, 



comparing SPS-mass estimates to total mass using dy- 
namics or lensing. Rather than exhaustively comparing 
all such results, we will focus on the studies that explic- 
itly derived S\mf for large samples of low- z ETGs. 

The s tudy most closely related to ours is iDutton et al.1 
(|2012aj) . who analyzed SDSS data for ETGs (colors and 
<t+, with SPS+ Jeans modeling). Their sample was larger 
but without the Sersic models and near-infrared photom- 
etry from SPIDER, which ought to provide more accu- 
rate results. Their results for a no-AC-NFW model are 
shown as a dotted orange curve in panel (a), which is 
reassuringly similar to our no-AC-NFW result, with the 
~ 30% residual difference in <5imf illustrating the level of 
systematic uncertainties for a fixed dataset and method. 
The solid line shows their refined result from a multi- 
parameter fit: they find DM halos that are slightly ex- 
panded, and consequently a slightly heavier IMF. 

In panel (b) we show results from the SAURON 
project, using both spectroscopically-based SPS models 
and detailed two-dimensional Jeans dynamical analyses. 
For the full ATL AS 3D sample, we used t he color-coding 
in figure 2(e) of iCappellari et al.l (|2012) to regenerate 
roughly their no-AC-NFW <5iMF-c e results. As a check, 
we a lso show their prototyp e individual galaxy results 
from Cappcllari ct al. (2006) (no-DM models), omitting 
the young galaxies that are thought to suffer from SPS 
modeling errors. As with SPIDER, the ATLAS 3D project 
found that the overall 5imf normalization was degener- 
ate to the DM assumptions, but the trend with a c was 
robust. Direct comparison in Figure 2] reveals that the 
two sets of SAURON results are consistent with SPIDER 
to within the errors, whi ch is also the case for a similar 
study of high-z galaxies (|Capperlari et al.ll2009f). 

P anel (c) show s resu lts from iThomas et all IpOlTh 
and iWegner et al.1 (|2012fl , who carried out spectroscopic 
SPS and detailed orbit-modeling of ETGs found in clus- 
ters (including Coma), using a variety of mass models 
(constant- Ydyn, no-AC-NFW, and cored halos). There 
is a systematic offset in 5imf between the two studies, 
but overall, the results are consistent with our no-AC- 
NFW results, both in amplitude and in slope. 

Panel (d) shows gravitational len sing results, prim arily 
from the SLACS survey of ETGs dTreu et alJl2010h . us- 
ing no-AC-NFW models as well as color-based SPS mod- 
els and stellar-dynamics constraints. Their <5imf normal- 
ization agrees well overall with ours, but their slope is 
somewhat steeper 

All four hybrid studies from the literature suggest the 
same basic result as ours: that the IMF of ETGs varies 
from MW-like at low a c to Salpeter-like at high cr e , mod- 
ulo some lingering uncertainties from the DM-IMF de- 
generacy. 

In addition to the large-sample studies, we exam- 
ine a few key single-galaxy results. These include 
extended kinematics data wi th NFW-based modeling 
(|Napolitano et al.l 120091 1201 If ) , as shown in panel (c) , 
and the "Jackpot" double lens (|Sonnenfeld et al.l l2012). 
in panel (d). These results are consistent with the trend 
from SPIDER. 

2 Many analyses of the SLACS data have been carried out in 
the literature, and not all of them are in agreement about the IMF 
conclusions. It is beyond the scope of this paper t o resolve these 
incons istencies, but we note that our own analysis in lTortora et al.1 
(2010) agrees with the SLACS-team analysis. 
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Figure 4. IMF mismatch-parameter vs. stellar velocity dispersion, comparing the SPIDER results (black curves) to various literature 
studies (see main text and legend, with abbreviations in the reference list). In several panels, typical error-bars are shown to the left. The 
filled CvD12 datapoint is from four galaxies' stacked spectra. There is generally good consistency in the inferred IMF trends from the 
different studies. 



Our culminating comparisons are with a completely 
different set of results, based purely on modeling of IMF- 
sensitive spectral lines (see Section rjj . Most of these 
lines are susceptible to degeneracies with elemental abun- 
dances (e.g., sodium or calcium), and we consider only 
those studies that have directly accounted for such ef- 
fects. 

We first show in panel (e) the results from CVD12, 
who studied a subset of the SAURON sample. They fit- 
ted spectral features across a wide wavelength range, fo- 
cusing on the IMF indicators Na I, Ca II, and the Wing- 
Ford band, while adopting a broken-power-law IMF form 
and fitting for the relevant elemental abundances. Their 
inferred <5imf values turn out to agree well with both 
the normalization and the trend versus a e from the SPI- 
DER results (for the no-AC-NFW models in particulars- 
Recalling that the high 5imf values from hybrid studies 
could be due to either a bottom-heavy (extra dwarfs) or 
a top-heavy (extra remnants) IMF, the CVD12 results 

3 CvD12 compared their T* results to total dynamical values 
from SAURON in order to check that they did not violate those 
constraints. However, they did not compare to decomposed dynam- 
ical inferences for consistency as we do here. 



agree with only the first solution. 

Note that the apertures probed here are different: R e /8 
for CvD12, and 0.3-0.7 R e for SPIDER (decreasing with 
<r*). The close agreement of the results on average thus 
implies that the IMF does not vary spatially on these 
scales, or that the AC model is the correct solution, and 
Simf decreases with galactocentric radius (which is plau- 
sible). 

Panel (f) shows results from lSmith et all (|2012f ). who 
studied a large sample of Coma-cluster ETGs. They used 
the same SPS models as CvD12 to fit near-infrared spec- 
troscopic line-indices, analyzing the Wing-Ford band 
and the Ca I line separately. The former line was suscep- 
tible to uncertainties in the Na abundance, but not the 
latter. We have converted their results to inferred <5imf 
by straightforward interpolation between the Chabrier, 
Salpeter, and a — 3 models in their figure 10. The final 
results are somewhat noisy and uncertain, but imply an 
overall Salpeter-like normalization, and no obvious trend 
with <t c . This agrees with all the aforementioned results 
at high-<7 e , but not at 1ow-<t c . It is possible that these 
Coma-cluster galaxies are genuinely different. We see 
very little dependence of <5imf on environment in SPI- 
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DER, but we also have very few cluster galaxies in our 
s ample. 

iSpiniello et~aT1 (|2012T > analyzed NaD, Na I, and Ti0 2 
lines in SDSS spectra of ETGs with > 200 km s _1 , 
comparing to the same SPS models as CvD12. They in- 
ferred a Salpeter-like IMF at low er*, which is consistent 
with SPIDER. At high er*, they inferred a ~ 3, which 
would imply Simf ~ 4 and violate the total mass con- 
straints both from SPIDER and from the lens galaxy 
that these authors also studied. This conflict suggests 
further work is needed on the line-index modeling. 

An ideal co mparison with our SP IDER results would 
be the work of iFerreras et al.l (|2012f ). who analyzed Na I 
line-strengths from the SPIDER parent dataset. Qualita- 
tively, their results imply a similar increase of <5imf with 
<r e , but more quantitative inferences will have to await 
more detailed analyses that allow for Na enhancement. 

It would also be fascinating to see whether or not 
spiral bulges follow t h e sam e <5tmf-ct* trends as ETGs. 
iFalcon-Barroso et al.1 ((2003) found that bulges showed 
anti-correlations between a* and Ca II line-strengths, 
similar to ETGs, but they could not determine if this 
was an IMF effect. We have attempted to investigate fur- 
ther, using (5tmf inference s for bul ges from the literatur e 
(e.g. Jde Blok et al.|[200l CvD12; lDutton et al.ll2012cl ). 
We find hints of an ETG-like Simf-v* correlation for the 
bulges, but we also find the same for the disks, which 
would be peculiar. The apparent scatter for both bulges 
and disks is enormous, suggesting more work is needed 
to address the systematic errors, and to understand any 
additional trends with detailed morphology. 

5. DISCUSSION AND CONCLUSIONS 

We have analyzed the dynamics and stellar populations 
of a large sample of ETGs from the SPIDER project, and 
found compelling evidence for heavier IMFs in the central 
regions of higher-ov galaxies. The IMF mismatch rela- 
tive to Chabrier is <5imf ~ 0.5-1.1 at a* ~ 125 km s _1 , 
and ~ 1.2-1.8 at at ~ 250 km s _1 . The <5imf~o* trend is 
robust to different modeling assumptions, and accounts 
for much of the tilt in the fundamental plane. The dis- 
tribution of DM is degenerate with the overall IMF nor- 
malization and difficult to constrain. 

We have performed the first general inventory of IMF 
results from a variety of studies in the literature, using 
both pure and hybrid techniques, and found that these 
generally agree well with our SPIDER results. There is 
remarkably widespread agreement on a Salpeter-like IMF 
for massive ETGs (cr c > 200 km s" 1 ). At lower a e , the 
data are still somewhat limited and the different results 
have not yet converged, but most of the studies point to 
MW-like IMFs. Agreement on a super-Salpeter IMF in 
the most massive galaxies (<j e > 275 km s -1 ) is also not 
yet universal. 

These results appear to be fully compatible with under- 
lying ACDM halos. However, more detailed conclusions 
about halo contraction or expansion are still elusive, and 
the data on galaxy centers do not clearly rule out alter- 
native DM models once the variable IMF is accounted 
for. 

More work is clearly needed to understand the system- 
atics in the different analyses; to build up better statistics 
on a wide range of galaxy types, environments, and red- 



shifts; and to determine which parameters correlate best 
with IMF va riations (e.g.. meta llicity or starburst inten- 
sity; CvD12; [Smith et alj|20l3 ). It' may also be partic- 
ularly helpful to venture beyond the centers of galaxies, 
using data from a wide baseline in radius to help break 
the IMF-DM degeneracies. 

It appears we are nearing convergence on determin- 
ing what the basic components of galaxies are (distribu- 
tions of stars and DM). The next challenge will be to un- 
derstand why these arrive at their distributions. What 
drives the power-spectrum in cloud fragmentation and 
star formation? How do baryonic processes interact with 
and re-shape their surrounding DM halos? 
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